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1 Introduction 

Motivated by the seminal theorem of Kharitonov on robust stability of interval polynomials H, 
a number of papers on robustness analysis of uncertain systems have been published in the past few 
years ||, ||, |[ 0, [j], ||, |[ [h]]. Kharitonov's theorem states that the Hurwitz stability of a real (or complex) 
interval polynomial family can be guaranteed by the Hurwitz stability of four (or eight) prescribed critical 
vertex polynomials in this family. This result is significant since it reduces checking stability of infinitely 
many polynomials to checking stability of finitely many polynomials, and the number of critical vertex 
polynomials need to be checked is independent of the order of the polynomial family. An important 
extension of Kharitonov's theorem is the edge theorem discovered by Bartlett, Hollot and Huang[Q. The 
edge theorem states that the stability of a polytope of polynomials can be guaranteed by the stability 
of its one-dimensional exposed edge polynomials. The significance of the edge theorem is that it allows 
some (affine) dependency among polynomial coefficients, and applies to more general stability regions, 
e.g., unit circle, left sector, shifted half plane, hyperbola region, etc. When the dependency among 
polynomial coefficients is nonlinear, however, Ackermann shows that checking a subset of a polynomial 
family generally can not guarantee the stability of the entire family jll], |l2|, [l3|| . 

Parallel to this line of research, robust stability of uncertain matrices has also received considerable 
attention. Bialas 'proved' that for robust Hurwitz stability of an interval matrix, it suffices to check all 
vertices . But Barmish and Hollot gave a counter-example to show that Bialas's claim is incorrect [l5| . 
Kokame and Mori considered Hurwitz stability of an interval polynomial matrix, and by using some result 
in signal processing theory, established a necessary and sufficient condition for robust stability fl^]. 

This paper studies robust D-stability of polytopic polynomial matrices, i.e., matrices with entries 
being polytopes of polynomials. We give a necessary and sufficient condition for robust Instability of 
Polytopic Polynomial Matrices, namely, the stability of a subset of this family guarantees the stability of 
the entire family. 



2 Preliminaries 



Definition 1 An interval polynomial matrix A — (pij) n xn is a matrix whose entries Py are interval 
polynomials, i.e., Pij — g° + qjjS + ... + <?™s m , q\- e k = 0, . . . , to, where k stands for 

superscript. A polytopic polynomial matrix A = (pij) n xn is a matrix whose entries py are polytopic 
polynomials, i.e., p VJ = YJk=i X ijPij> X ij > °> EfeLi X ij = M = 1, where p^ are fixed 

polynomials. 
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Definition 2 Suppose fl C R" +1 is an m-dimensional polytope. Its supporting plane H is defined 
as an n-dimcnsional affine set, satisfying Q n H =/= 0, and all points of Q lie on the same side of H: Its 
exposed set is denned as the intersection of f2 and its supporting plane H; Its exposed edge set is defined 
as the one-dimensional exposed set. 

An n-th order polynomial can be regarded as a point in its (n + l)-dimensional coefficient space. 

Definition 3 Given an open region D in the complex plane, the polynomial matrix A is said to be 
Instable, if all roots of detA = lie within D; A polynomial matrix set A is said to be D-stable, if every 
member in A is D-stable. When D is taken as the open left half of the complex plane, Instable is also 
called Hurwitz stable. 

In fact, Edge Theorem holds for more general stability regions. For simplicity, we only consider 
simply-connected stability regions in this paper. 

Definition 4 Given an interval polynomial set J-~(s) — {5Z™ li^i 1i e [<7-i9i]}, its Kharitonov 
vertex set is K% = {ft, /f, /* /£}, and E% - {A/| + (1 - A)/*, (s, t) G {(1, 2), (2, 4), (4, 3), (3, 1)}, A G 
[0, 1]} is called its Kharitonov exposed edge set, where 

fk =%+ lx 3 + ?2« 2 + • • • fk=2o + 9i* + ?2« 2 + ■ ■ ■ 
fk =%+ + g 2 s2 + ■■■ fk=% + «i* + g 2 s2 + ■■■ 

Consider the polytopic polynomial sets 

■Pij = {ELi ■ A* > 0, ELi 4 = 1} (m = 1, . . . , n) 

p%j are fixed polynomials, k = 1, . . . , m. 

Their vertex sets are 

Kij = {Ptj k = l,...,m} {i, j = 1, . . . , n) 
and by definition, their exposed edge sets are contained in 

Eij = {Xp-j + (1 - X)Pij, s,t = l,...,m} = 1, . . . ,n) 

Let A = {(Pij) nxn '■ Pij S Pij > j ' — 1? ■ • • i n }? and let P™ be the set of all permutations of 1, 2, . . . , n. 
Definition 5 Define e_4 as 

//„ \ . Psi s EE s i s ,(h,...,l n ) eP™,s = l,...,n \ , , 

<, {Pijjnxn G K ■ i = 1 Z - 1 Z + 1 n I ^ J 

It is easy to see that, is produced by taking only one entry from its exposed edge set in every 
row/column and all other entries from their vertex sets in A. 



3 Main Results 

3.1 Polytopic Polynomial Matrices 

Consider the polytopic polynomial matrix set 

-4 = {(Pij)nxn ■ where p y - G Vij, i,j = 1, . . . ,n} (3) 

Suppose VA G A, deg(detA) — const. 

Theorem 1 ^4 is D-stable if and only if is Instable. 

Proof: Necessity is obvious. To prove sufficiency, suppose ej± is D-stable. Let 
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By using Laplace formula on the first column, we have 

detA = piiMn + . . . + p n iM nl 

Let 



Pn P12 



T 



Pal Pn2 ■ 

It is easy to see that A G T and VTeT 



pi„ \ and p* a e Vn ; i = 1, . . . , n 

\ , Pij are entries of A 
Pnn ) j = 2, . . . , n 



detT = p* u M n + ■■■+ P* nl M nl 
Apparently, detT is an affine function of . . . ,Pn\- By Edge Theorem 

detT is D-stable <^> the edge set of detT is D-stable. 

The edge set of detT is 

EuMa + K n M n, i = 1, ■ • ■ ,n ^ . 



The corresponding matrix collection is 
A\ = < 



qil (Pst)nx(n-l) 



fti G ^i;« G {1, . ..,n} 
9fei G -ftTfci;* = 2, . . . ,n 
Pfet ; Pit are entries of A 
k = 1, . . . ,i — 1, i + 1, . . . , n 



In this case, T is instable <^> .Ai is D-stable. Moreover, V-Ai £ ,4i, there exists i 6 {1, . . . ,n},qu E 
En, qui G iffei, k = 1, . . . ,i — 1, i + 1, . . . , n such that 



/ 9n P12 



A x = 



Qil Pi2 (Pst)nx(n-2) 
\ q n l Pn2 

Again, by using Laplace formula on the second column, we have 

detAi = P12M12 + . . . + Pn2M n2 

Set 

911 P\2 Pl3 •■■ Pin 
Qnl Pn2 Pn3 ■■■ Pnn 

Let T\ = lJ Aig ^ 4i B, then, its edge set A2 is 



B= I 



qij are entries of A\ 
Pij are entries of A 

Pij G Vij 
i,j = l,...,n 



' / 911 <7i2 
■ • • q 3 2 

\ (Zral 9n2 



(Pat) 



\ 3 G {l,...,n};/ = l,...,n 
m = 1,1 ^ i 
m = 2,l^j 



K, 



Ira 1 



qim g 



m = 1,1 = i 



m = 2,l = j 
Pit are entries of A 
t = 3, . . . , n 
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and A\ C T\. By definition and Edge Theorem 

71 is Instable Ai is D-stable. 
By repeating the process above, we have 

A2 , • ■ ■ 5 An—1 , A n ', T~2, ■ ■ ■ , %i-l 

Tk is D-stablc <^> Ak+i is D-stable 

Ak is D-stablc <= Tk is D-stable k = 2, . . . , n — 1 



where Ak is the collection of 



< 



qi k i 

k V 9nl 



lilk 

■ ■ ■ (Pst) 
<li k k 

Ink 



qit e 



Ei t t,l — *t 

K itt ,l^ i t 
ii,...,i k e {l,...,n} 
t= l,...,k 
I = 1, . . . , n 
Pis are entries of A 
s = k + 1 , . . . , n 



Thus, for each element of Ak, its entries have the following characteristics: for the first k columns, all 
entries of each column belong to their vertex sets except that one entry belongs to its exposed edge set, 
and the entries of the remaining n — k columns are the corresponding entries of A. Hence, VA n e A n , we 
have 

'/ qn ■■■ qin \ 



A n < 



Qhi 



qi n n 



qit e 



E itt ,l = it 
K itt ,l # i t 
h,.. .,i n G {1,. ..,n} 
i = 1, . . . , n 
i = 1, . . . , n 



If « s = i 4 for some pair i s , i t , without loss of generality, suppose i\ = 12 = 1, namely 

qn qu ■■■ qin 



A n — 



, 9n G En, qi2 G E 



12 



<7nl <Zn2 • • ■ q-nn 

By using Laplace formula on the first row of A n , we have 



detA, n = quMn + q 12 M 12 + ^ q u M u 



i=3 



By Edge Theorem 

A n is instable ^ quM n + ql 2 M 12 + £" =3 <7h m h and 

n 

g^Mn + (712M12 + ^2 QuM u are D-stable. 
The corresponding matrices are 

/ 0?n O19 ... Ol„ \ 

,<?ii G K 1U 
,Qi2 G #12, 



qu 


qi2 ■ 


■ qin 


q n i 


q n 2 ■ 




qu 


q\2 ■ 


■ qin 


q n i 


qn2 ■ 





4 



which belong to 



< 



Qhi 



{ V Qui 

So, An is D-stablc e_4 is D-stablc 



Qln \ 

qi in 



Qnn 

Thus 



qi sS e £ 4sS , (ii, • • .,«„) e f£ 

Z = 1, . . . , n 
s = 1, . . . , n 



> . 



e^4 is D-stable 44> .A„ is £> — stable 44> 7^_i is D-stablc 

=>- -4n_i is D-stablc . . . % is D-stable 

=> Ai is D-stable <^> T is D-stablc =>■ A is D-stable. 

That is to say, ej, is D-stablc =>■ \/A e .4, A is D-stable, namely, A is D-stable. This completes the proof. 
Remark 1 When m = 2, i.e. A = (Pij) nxn , where 



Pij = Pi, + Kp\ 3 K g [o, i] 

Pij ) are fixed polynomials i , j 



1 ,n 



Let 



Bi = {(p. 



Pij satisfies (4)} 



(4) 



(5) 



In this case, the vertex sets of p^ are , _p° 3 } , and their exposed sets are exactly themselves, namely, 
{p° + Xijpjj}. The corresponding conclusion has more concise form, this is due to the simplification of 
the edge set of B\ 

(Zi,...,Z„) GP£;s = l,...,n 

e Sl = \ (Pij)nxn x / [0, 1] i = l s 

" V < {0,1} i + l a 



3.2 Interval Polynomial Matrices 

Consider the subset of R nxn (s) 

&2 = {(Pij)nxn, Pij are interval polynomials } 

where R nxn (s) is the collection ofnxn polynomial matrices. Assume eg 2 is Hurwitz stable. 
Consider detA, using Laplace formula on first column, we have 



(6) 



detA = puM n 



■PnlM„i 



Similar to the proof of theorem 1, and by resort to the Generalized Kharitonov Theorem, we have 

Theorem 2 If VA 6 B2,deg(detA) — m. Then £> 2 is Hurwitz stable eg 2 is Hurwitz stable, 
where eg 2 i s 

Pih G E% > ;i = 1,. • • ,n 

(Pij)nxn ■ Psi„ G K s i s ', l s ^ i s = 1, . . . , W, 

(Zi, ...,l„)e^;s = l,...,n 
and , are defined in Definition 4. 

Remark 2 Theorems 1 and 2 can be viewed as a generalization of the Edge Theorem and Kharitonov 
Theorem to MIMO case. Stability test of the entire family is reduced to a critical low-dimensional subset. 
No extra lemma from signal processing is needed in our proof. Furthermore, our results can be easily 
extended to polynomial matrices with complex coefficients. 
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